Theory of optical spectral weights in Mott insulators with orbital degrees of freedom 
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Introducing partial sum rules for the optical multiplet transitions, we outline a unified approach 
to magnetic and optical properties of strongly correlated transition metal oxides. On the example of 
LaVC>3 we demonstrate how the temperature and polarization dependences of different components 
of the optical multiplet are determined by the underlying spin and orbital correlations dictated by 
the low-energy superexchange Hamiltonian. Thereby the optical data provides deep insight into the 
complex spin-orbital physics and the role played by orbital fluctuations. 

PACS numbers: 75.30.Et, 75.10.-b, 75.10.Jm, 78.20.-e 



I. INTRODUCTION 

Charge localization in Mott insulators is not perfect 
as electrons still undergo virtual transitions to neigh- 
boring sites in order to retain partially their kinetic en- 
ergy. These high-energy virtual transitions across the 
Mott-Hubbard gap ~ U are crucial for magnetism - 
this quantum charge motion leads to superexchange (SE) 
interactions* between local degrees of freedom. It is fre- 
quently not realized, however, that the same charge ex- 
citations are responsible in Mott insulators for the low- 
energy optical absorption. Therefore, the intensity of the 
optical absorption and the SE energy are intimately re- 
lated to each other via the well known optical sum rule^ 
which links the integrated optical conductivity and the 
kinetic energy. In Mott insulators the latter is due to vir- 
tual exchange processes and hence the thermal evolution 
of spectral weights and SE energy are related^ 

A qualitatively new situation is encountered in Mott 
insulators with partly filled d orbitals. Because of or- 
bital degeneracy, virtual charge excitations that are seen 
in optics reflect the rich multiplet structure of a tran- 
sition metal ion determined by Hund's exchange cou- 
pling Jjj, and orbital degrees of freedom contribute then 
to the SEi4 When spin and orbital correlations change 
the individual components of the optical multiplet re- 
flect characteristic spectral weight transfer. Moreover, 
the cubic symmetry is spontaneously broken by orbital 
and spin order, and thus one expects anisotropic opti- 
cal absorption. Indeed, pronounced anisotropy was re- 
ported for LaMn03;£> both for the A- type antiferromag- 
netic (AF) phase 6 and for the orbital ordered phase above 
the Neel temperature TV- Recently, the anisotropy in 
optical absorption and its strong temperature depen- 
dence near the magnetic transitions were found for cu- 
bic vanadates^ This latter example is even more puz- 



zling as the magnetic properties are anomalous^ and 
neutron scattering experiments^ have revealed nontriv- 
ial quasi one-dimensional (ID) correlations of spin and 
orbital degrees of freedom that are surprising for crystals 
with nearly cubic symmetry. Indeed, a theory of spin 
and orbital states in cubic vanadates predicted quasi ID 
spin-orbital correlations due to a spontaneous breaking 
of the cubic symmetry in the SE models 

It is our aim to outline a unified picture that links opti- 
cal and magnetic properties at orbital degeneracy. Start- 
ing from the low-energy spin-orbital model we derive par- 
tial sum rules for the different excited states. Thereby 
we provide a rigorous theoretical basis for the analysis of 
optical spectral weights and show how the evolution of 
magnetic coherence manifests itself in optics as intensity 
transfer between different excitations (upper Hubbard 
bands). This explains the origin of dramatic variation of 
the optical absorption and its anisotropy with temper- 
ature T in manganites^ vanadates^ and ruthenates^i 
where at low T the high-spin band carries most spec- 
tral weight for the directions with ferromagnetic (FM) 
spin correlations. We illustrate this idea for the case of 
C-type AF (C-AF) phase 6 - of cubic vanadates with de- 
generate tig orbitals. We show that the predicted quasi 
ID spin-orbital correlations r°« realized in C-AF phase 
of LaV03, are reflected in the T-dependence of optical 
weights derived from the SE model. 

The paper is organized as follows. First we present a 
generic structure of the SE interactions in a Mott insula- 
tor with orbital degrees of freedom in Sec. II, and relate 
them to the intensities in optical absorption. On the ex- 
ample of the SE interactions encountered in LaV03, we 
analyze next spin, orbital, and joint spin-and-orbital cor- 
relations which determine the optical intensities. In this 
way, we arrive at a set of self-consistent equations which 
are solved in Sec. Ill, where we present the numerical 
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results for the optical sum rules for individual high-spin 
and low-spin excitations, as well as for the spin and or- 
bital SE interactions. These results and their comparison 
with experiment are discussed in Sec. IV, where we also 
give a summary and more general conclusions. 



II. THEORY 

The SE interaction in a cubic Mott insulator with ions 
having orbital degrees of freedom has a generic form, 



Hj = H 8 



<*i>llT 



and consists of separate spin (H s ) and orbital (H T ) in- 
teractions, and of a dynamical coupling between them 
(H ST ). This complex form of Hj, given by Eq. QJ, fol- 
lows from the terms H^\ij) for each bond (ij) along a 
given cubic axis 7 = a, b, c, arising from the transitions 
to various upper Hubbard bands labelled by n. The opti- 
cal intensity of each band n, for the photon polarization 
along a cubic axis 7, is determined by the respective SE 
energy: 

t2 /.oc . 

"" / a^( U )d U = -^ = - V (^m. (2) 



Here, ao is the distance between magnetic ions (the tight- 
binding model is implied), and (^Hn (ij)) is the SE in- 
teraction for a bond (ij) along axis 7. The first equality 
in Eq. J5J follows from the optical sum rule for a given 
transition n, and relates the kinetic energy to the 
optical conductivity an(w) for this band, while the sec- 
ond equality relates the associated kinetic energy to the 
SE energy via the Hellman-Feynman theorem^ 

Experimental data is often presented in terms of an 
effective carrier number [see, e.g., Eq. (2) of Ref. 0, 

^eSn = (2™o«o/ 7r e 2 ) Jq 00 Un\uj)duj, where m is the free 
electron mass, and vq = af, is the volume per magnetic 
ion. This gives an optical sum rule as follows: 

Each level n of the multiplet represents an upper Hub- 
bard band with its own spin and orbital quantum num- 
bers. The key point is that while full kinetic energy and 
corresponding total intensity may show only modest In- 
dependence and almost no anisotropy, the behavior of the 
individual transitions is much richer, and directly reflects 
the ground state correlations via the spin and orbital se- 
lection rules. 

Hund's exchange separates the lowest (high-spin) ex- 
citation from the next (low-spin) one by: ~ 5 Jr 
in manganites, 12 3J# in vanadates, 10 and 2Jjj in 
titanates. 1 - When the high-spin excitation, broadened by 
its propagation in crystal and by many-body effects, has 



a smaller linewidth than the above multiplet splitting, 
it may show up in optical spectroscopy as a separate 
band. This is in fact nearly satisfied for typical values 
of J H ~ 0.6 - 0.7 eV and hoppings t ~ 0.4 eV (~ 0.2 eV) 
for e g (t2g) orbitals. 14 Higher bands overlap and mix up 
with d — p transitions, though. 

In a particular case of vanadates, one has three optical 
bands n = 1,2,3 arising from the transitions to: (i) a 
high-spin state 4 A 2 at energy U—3Jh, (ii) two degenerate 
low-spin states 2 Ti and 2 E at U, and (hi) 2 T 2 low-spin 
state at U + 2Jh*^ Using 77 = Jh/U we parametrize 
this multiplet structure by: R = 1/(1 — 3rj) and r — 
1/(1 + 2rj). In LaVC>3 xy orbitals are singly occupied, 8 
and one obtains a high-spin contribution H{ c \ij) for a 
bond (ij) along c axis: 

h[ c) = —JRfSi-Sj +2) (| - n-rj), (4) 

while for a bond in an (a, b) plane: 
1 
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(ab) 



-Jr(s v S 3 +2)(\-t*t*). 



(5) 



In Eq. Q pseudospin operators fi describe low-energy 
dynamics of (initially degenerate) xz and yz orbital dou- 
blet at site i; this dynamics is quenched in i?j afc ' 1 JSJ. Here 
■g(Si ■ Sj + 2) is the projection operator on the high-spin 

state for S = 1 spins. The terms Hn\ij) for low-spin 
excitations (n = 2, 3) contain instead the spin operator 
(1 — Si ■ Sj) (which guarantees that these terms vanish for 
fully polarized spins on a considered bond, (Si ■ Sj) = 1): 
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while again the terms H^ b \ij) differ from Hn\ij) only 
by orbital operators: 
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where upper (lower) sign corresponds to a(6)-axis bonds. 

First we present a mean-field (MF) approximation for 
the spin and orbital bond correlations which are deter- 
mined self-consistently after decoupling them from each 
other in Tij (Q. Spin interactions, 

Hs = Jab E Si ' S 3 ~ J c E Si ' ( 8 ) 

(ij)ab (ij)c 

depend on exchange constants: 
1 



J' 



J, 



ah 



r]R- {R-i 1 R- W ){\ + (r l -T 1 ))-2r 1 r(T?Tj) 
l-r)R-r)r+ (R-r)R-r]r)(l + (rfr/))], (9) 
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determined by orbital correlations. In the orbital sector 
one finds 

H r = E ^ • f i J ^ ^rrfr]] + J T ah £ rfr/, 

(10) 

with: 



2 

JT I J 

J ab — A J 



(l + s c )R+(l-s c )r)(R + r) 
(l-s ab )R+(l + s ab )r](R + r) , (11) 



depending on spin correlations: s c = (Sj • Sj) c and s a fc = 
— (Si ■ Sj) a b- In a classical C-AF state (s c = s ab = 1), 
this MF procedure becomes exact, and the orbital prob- 
lem maps to Heisenberg pseudospin chains along c axis, 
weakly coupled (as 77 <§; 1) along a and b bonds, 



Hi® = JR 
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(12) 



releasing large zero-point energy. Thus, spin C-AF order 
and quasi ID quantum orbital fluctuations support each 
otherifi 

In addition to spin-orbital SE Hj JT} , orbitally degen- 
erate systems experience the Jahn- Teller (JT) interac- 
tions — coupling of orbitals to lattice distortions may 
lead to a structural transition, lifting orbital degeneracy 
by the term, 
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Ti T., 



(13) 



in the present vanadate model. 10 Interactions V a b > 
originate from the coupling of nearest-neighbor t 2g or- 
bitals in (a, b) planes to the bond stretching oxygen vi- 
brations in corner-shared perovskite structure. They gen- 
erate antidistortive oxygen displacements and staggered 
orbital order (supporting SE), whereas the V c > term 
due to the GdFeC>3-type distortion^ favors ferro-orbital 
alignment along c axis, and thus competes with SE. 

The complete model H = Ttj + TLv represents a non- 
trivial many-body problem. Interactions are highly frus- 
trated, leading to strong competition between different 
spin and orbital states. We leave this complex problem 
for a future study, and present here the conceptually sim- 
pler case of LaVC>3 with C-AF order. In this phase spins 
are FM along the c axis at low T, and orbitals fluctuate 
on their own. This justifies a posteriori a perturbative 
treatment of joint spin-orbital correlations and allows one 
to determine them in a simple analytical way. 

We begin with the orbital Hamiltonian H T + Hv [Eqs. 
l[TU)l and (^21] which has the form of an XYZ model. 
As interchain couplings are weak and of t'tj form, the 
problem is best handled by employing Jordan- Wigner 
fermion representation^ After decoupling t? t? terms in 
fermionic density and bond-order terms, one finds that 



the staggered orbital order parameter r = |(t*)|, the or- 
bital ordering temperature T T , and the temperature de- 
pendence of orbital correlations: (t* Tj') c = — 

( T i T j)c = — K 2 , (TiTj)ab = ~T 2 , do follow from twO 

self-consistent equations: 



E 



\2s k J 
( Jc 



tanh ( — ) . 
\2TJ 



(14) 



« = E(^)^tanh@), (15) 

k 



where = [J 2 cos 2 k + (h T ) 2 ] 1 / 2 is the ID orbiton dis- 
persion, J c = J T C + 2k( Jl - Vc), and h T = 2r( J C T + 2 J T ah + 
2V a b — V c ) is the effective field. We set fc# = 1. 

The short -range spin correlations s-y determine J^ 
(|ll|l and are finite also above Tn- We derived them 
by solving exactly a single bond (ij) within the mean- 
field cx (S z )r^ originating from neighboring spins. This 
is the simplest cluster mean-field theory known in the 
theory of magnetism as Oguchi method^ For a FM 
bond along c axis one finds even an analytic solution: 
s c = (Z -Z 1 -2Z 2 )/Z, where Z = 1+2 coshx+2 cosh 2z, 
Z x = (1 + 2 cosh x) exp(-2 J*/T), Z 2 = exp(-3J*/T), 
with Z = Z + Z 1 + Z 2 , x = h s /T, h s = (J s c +U° b )(S z ). 
The s a b correlation function for an AF bond can be found 
numerically. 

Now we turn to the dynamical coupling between orbital 
and spin sectors, denoted as H ST term in Eq. (I}. In the 
present case of C-AF ground state it contains mainly 
contributions due to c axis bonds, and reads as follows: 



*E 

to% 



6 (Si ■ Sj)S(fi 



(16) 



Here, K = \ j(R - r)R - r/r) and 6(A) = A — (A) im- 
plies the fluctuating part of an operator, which goes be- 
yond the MF decoupling. Treating H ST within the high- 
temperature expansion, we found that the joint spin-and- 
orbital correlations, = (6 (Si ■ Sj)5(fi ■ Tj)), are given 



as follows: 



16T (((^ 



s,) 2 ) - (Si ■ SjY 



(17) 



with ((Si ■ Sj) 2 ) = 1 + 3Z 2 /Z in the present Oguchi ap- 
proximation. This high-temperature expansion for /y is 
valid when spin fluctuations are weak as in the C-AF 
phase, and one finds that the joint correlations van- 
ish at T — > when fully polarized spins decouple from 
the orbital sector. 



III. NUMERICAL RESULTS 

Taking the SE energy scale J ~ 40 meVjS we set 
Tyv = 0.4J, while r\ ~ 0.12 follows from spectroscopy^ 
It is quite natural in t 2g systems that order and disorder 
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FIG. 



1: (color online) (a) Kinetic energy K[ c ^ for the high- 
spin excitations within classical and quantum models. Solid 
(dashed) line with (without) joint spin-orbital fluctuations, 
(b) Intersite correlations along c axis: spin s c = (Si-Sj) c , or- 
bital (fi ■ Tj) c , and spin-orbital fij. Orbital order parameter 
r is also shown. Dotted line (filled circles) for —fij obtained 
from high-temperature expansion 117H (by exact diagonaliza- 
tion). Parameters: rj = 0.12, V c = 0.9J, V ab = O.2.J. 



in the spin and orbital sectors support each other; 10 ! 19 
and indeed T T ~ T N in LaV0 3 i While the microscopic 
reasons of this behavior are subtle, we control T T by tun- 
ing lattice mediated couplings Vy — we found that the 
orbital transition occurs near TV for: V c — 0.9J and 
V a b = 0.2J. Then the energy of the G-AF phase is still 
lower than that of the G-AF phase, but a slight increase 
of Vy (by a factor 1.3) due to a stronger tilting of VC>6 
octahedra gives the G-AF phase, observed in YVO31& 

Following Eqs. @, {TTJ, lfTl)l. (|T3|l . and O, we cal- 
culated orbital, spin, and joint spin-and-orbital correla- 
tions, and next used them to determine the kinetic en- 
ergies Kn^ © due to each Hubbard subband. For in- 
stance, 



K[ C) /2J : 



1 



-I? 



Si • Sj 



h 



(18) 



Quantum effects beyond the MF theory are particularly 
pronounced in 
approach —K± 



(c) 



■K[ c ' [see Fig. ^a)]. In the classical 
l|18fl increases with decreasing T only 
if T < T N (when both r > and (S z ) > 0). This is 
qualitatively different in the quantum model when the 
orbitals fluctuate, the orbital order parameter r is not 
more than half of its classical value, and the spin correla- 
tions s c are finite above TV and decay slowly for T > Tn 
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FIG. 2: (color online) Kinetic energy K}?' (solid lines) and 
total (dashed lines) in units of 2J for: (a) c axis and 

(b) afe-plane polarization. Parameters as in Fig. 1. Filled 
squares in panel (a) represent the effective carrier number 
in LaVOa which includes the sum of the peaks 1 and 2 
below 3 eV in Fig. 3 of Ref.0. The experimental data follows 
well the calculated intensity of the high-spin transition, n = 1. 



[Fig. Gib)]. In this case — K[ c ' is enhanced at T = 
by the orbital fluctuations (fi ■ fj) c ~ —0.43, being close 
to those found for a ID AF spin 1/2 chain. Kinetic en- 
ergy gain [Eq. I|18l) ] gradually decreases with increasing 
T, and is reduced by half at T ~ 2T/v from its value 
at T = [Fig. nja)]. We also note that joint spin-and- 
orbital fluctuations /y develop at finite T and contribute 
significantly at T > T^iSi An opposite behavior occurs 
in the limit of large JT orbital splitting that freezes out 
orbital fluctuations^ — then the temperature variation 
and the anisotropy of optical absorption are quenched, 
and (fi ■ fj)y approaches its classical limit (—0.25) in all 
directions. 

' (7) @ exhibit pro- 



The optical intensities 



"cff.i 



nounced anisotropy between c and ab polarizations (Fig. 
|2J), particularly those of high-spin transitions (n = 1) 
at low energy U — 3Jh- These low-energy intensities be- 
have here in opposite way for c and ab polarizations when 
temperature increases and the spectral weight is trans- 
ferred between the high-spin and low-spin bands. The 
total intensities have a much weaker temperature depen- 
dence than individual contributions. The theory repro- 
duces quite well the observed^ variation of N^g with T 
[Fig. Ufa)]. We recall that the temperature variation 

(c) (c) 

of the kinetic energy K{ ' and N^J is due to evolution 
of both spin and orbital correlations; spin-only ordering 
with frozen orbitals cannot explain a factor of two en- 
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FIG. 3: (color online) Exchange constants as functions of T 
along c (a/6) axis: (a) orbital Jj ill 1 1 . and (b) spin 0. 
Parameters: 77 = 0.12, 14 = 0.9 J, V ab = 0.2 J. 



hancement of iV Gfr in a range of T < 2Tjy. Therefore, 
the optical data of Ref. indicate that the JT orbital 
splitting in LaVC>3 cannot be largeiSi 

The microscopic reasons of anisotropy in the optical 
absorption are revealed by studying the effective ex- 
change constants, given by Eqs. @ and ijTT|l . While 
the spin sector is always 3D, the orbital one shows a di- 
mensional crossover from 3D to quasi ID correlations 
when the C-AF order develops [Fig. 02a)]. The orbital 
singlet correlations along the c axis enhance strongly FM 
J* at low T. Unlike in pure spin systems, the exchange 
interactions are temperature dependent [Fig. 13(b)], 
and J* decreases fast with decreasing intersite orbital 
correlations [Fig. IHb)]. Only at T ^> Tjv it is much 
weaker than the AF one, in agreement with conventional 
Goodenough-Kanamori rules. 



IV. DISCUSSION AND SUMMARY 



The basic experimental findings in the optical spectra 
of LaVC>3^ such as: (i) pronounced temperature depen- 
dence of c axis intensity (changing by a factor of two 
below 300 K), (ii) large anisotropy between the optical 
spectral weights along c and a/b axis (both below and 
above T)v), are qualitatively reproduced by our theory 
(see Fig. El ^ This strongly supports the picture of quan- 
tum orbital chains in the C-AF phase of vanadates. 

The spectral shape of optical absorption in LaVC>3 is 
highly intriguing, developing sharp coherent peak be- 
low 3 eV at T < Tn 7 < When its two-peak structure 



is interpreted^ as following from the multiplet splitting 
~ 3Jff, the value of Jh would be seriously underesti- 
mated, giving Jh — 0.21 eV, which is much smaller that 
the respective value (Kanamori parameter) 0.68 eV in 
the free ioniii Therefore, we suggest that the absorption 
below 3 eV is arising entirely from high-spin band cen- 
tered at U — 3Jh- The coherent peak emerging at low 
temperature could be then interpreted as a bound state, 
similar to what occurs in ID systems^ This is quite nat- 
ural when spins along the c axis are fully polarized and 
coherent orbital chains are formed. 

While the precise structure of the optical band requires 
further work, we may already obtain some useful infor- 
mation. Taking the width of 8t for the optical band^ 
and fitting the position and width of the spectral den- 
sity in Fig. 2(a) of Ref. Q, one finds U - 3J H ~ 2.3 eV, 
t ~ 0.2 eV, and thus JR ~ 70 meV. As higher bands with 
n = 2, 3 are not resolved in the experiment, we cannot 
fix Jh directly from the multiplet structure, but consid- 
ering r\ = Jh /U ~ 0.13 we reproduce basic energy scale 
J ~ 42 meV, being close to J ~ 40 meV deduced from 
neutron scattering data»^ In fact, we can also determine 
J directly from the comparison of the theoretical result 



-K[ c) ~ 1.1(2J), and the observed N^> ~ 0.21 at T = 
[see Fig. 0a)]. Using Eq. © with a Q = 3.91 Afi one 
obtains then J ~ 48 meV, a value remarkably consistent 
with both above estimates. Finally, we remark that the 
present interpretation of the experimental data of Ref. 
gives (at n ~ 0.13) Jh — 0.5 eV and U ~ 3.8 eV, which 
are somewhat lower than the atomic values^* This reduc- 
tion may be attributed to covalency and/or many-body 
screening effects in a solid ^ 

Summarizing, we proposed a new approach employing 
partial optical sum rules in Mott insulators with orbital 
degeneracy, which provides a theoretical framework for 
common understanding of the optical and magnetic ex- 
periments, both determined by the superexchange. Con- 
sidering the example of LaV03 with C-AF order, we 
have shown that pronounced temperature dependence 
and strong anisotropy of the optical absorption indeed 
follow from quasi ID quantum spin-orbital correlations, 
being radically different from the classical expectations. 
A satisfactory agreement between the values of J ex- 
tracted independently from the magnetic and optical 
data in LaV03 demonstrates that superexchange inter- 
actions are indeed responsible for the distribution of the 
optical spectral weight in Mott insulators. 
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